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Abstract 

The scalar two-loop master diagram is revisited in the massive cases needed for the 
computation of boson and fermion propagators in QED and QCD. By means of the 
causal method it is possible in a straightforward manner to express the propagators as 
double integrals. In the case of vacuum polarization both integrations can be carried 
out in terms of polylogarithms, whereas the last integral in the fermion propagator 
cannot be expressed by known special functions. The advantage of the method in 
comparison with Feynman integral calculations is indicated. 
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1 Introduction 



It is the big advantage of Epstein and Glaser's causal approach to quantum field theory [1] 
that it explicitly uses the causal structure of the theory, which is invisible in the traditional 
Lagrangian approach using Feynman rules or in path-integral methods. As a consequence, 
it is possible to prove the general properties of the S-matrix in very direct and transparent 
manner [2]. However, the full strength of the causal method comes out in higher order 
calculations and can be summarized in the following points: 

1. ) There exists no ultraviolet problem in the causal theory, so that no regular ization 
is necessary. 

2. ) All integrals are four-dimensional, so that no 7^-problem exists. Furthermore, 
the number of non-trivial integrations is reduced to the minimum, essentially one one- 
dimensional integral for every additional loop. 

3. ) The method is inductive in the order n of perturbation theory, so that the work 
done in lower orders is completely and automatically utilized in higher orders. 

The main reason why the Feynman rules are not optimal for loop diagrams is the 
following. Suppose that time ordering can simply be done by multiplication with step 
functions @{x^ — x^), as in the scalar theory we are going to consider. Then, in calculating 
a time-ordered product according to the rules, every pairing of field operators gets such 
a 0-function which leads to the Feynman propagators Dpixk — xj). However, the time 
ordering of the vertices is already specified by less G-functions. Without the temporal 
0-function the pairing function is instead of Dp. Consequently, some propagators 
Dp can actually be simplified into D^ which is simpler to integrate because it contains a 
(5-distribution in p-space. For example, in the causal theory the third order scalar vertex 
function is given by 

A(xi,X2,X3) = Df{xi - xs)D+(x3 - X2)Df{xi - X2) - DpD^^'D- + D+D"-''D^^K (1.1) 

If one substitutes D^ = Dp — D"''" in the first term, one arrives at the usual Feynman 
form 

A = DpDpDp - DpD'''"{Dp + D') + D+D^'^'D^^^ (1.2) 

because the last two terms are equal to D"''" D^'" D^^^ which vanishes by the support prop- 
erties of the advanced and retarded distributions. But (1.1) is simpler to calculate than 
(1.2) because every member contains one D^. This advantage is strongly increasing in 
higher loop diagrams. 

It is the purpose of this paper to illustrate these features for the so-called two-loop 
master diagram of Fig.l. Interestingly enough, this diagram was already computed with 
the causal method by various authors without knowing it. The first were Kallen and Sabry 
[3], their work was extended by Broadhurst [4] to other mass cases. The most extensive 
higher order calculations using the " dispersive method" were carried out over many years 
by the Italian group Mignaco, Remiddi [5], Barbieri [6] and others (see [7] and references 
given therein). All these authors base their calculations on analytic properties of Feynman 
integrals which are often referred to as Cutcosky rules. The lack of understanding of these 
analytic properties has created the problem of anomalous thresholds. Fortunately, the 
latter do not appear in diagrams with two and three external legs, but they do appear in 
four- and more legs-diagrams [8]. 
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Such basic subjects are discussed in the following section. As a first application, we 
then briefly describe the calculation of the scalar vertex function with arbitrary masses 
which is needed at various places in the later two-loop calculations. In Sect. 3 we turn 
to the master diagram for vacuum polarization. We show the calculation of the causal 
distribution in some detail to localize the infrared divergences which appear in the case 
of vanishing photon mass. The most difficult integration is the dispersion integral for the 
splitting. The necessary techniques to handle polylogarithms are shown in Sect. 4. Finally 
we consider the more complicated case of the fermion propagator. By ingenious tricks, 
Broadhurst [4] succeeded in expressing this propagator by a one-dimensional integral over 
the complete elliptic integral of the first kind K(k) together with elementary functions 
which can easily be computed numerically. In the causal theory one needs not be a 
genius. We get the result by a slight change of the standard procedure. As far as the spin 
structure is concerned, it is known that the general two-legs diagram can be expressed in 
terms of the scalar two-point functions [9]. Therefore, we restrict ourselves to the scalar 
case here. 



2 The Causal Method and the Scalar Vertex 



We first summarize the main ingredients of the causal method, for details we refer to 
[1,2]. In the causal theory the S-matrix is viewed as an operator- valued distribution of the 
following form 

S'(5) = 1 + X!~F dxi...dxnTn{xi,...Xn)g{xi) ■ ...g{xn), (2.1) 

-, Tl> J 

n=l 

where g € S'(R^), the Schwartz space of functions of rapid decrease. The test function g 
plays the role of " adiabatic switching" and provides a cutoff in the long-range part of the 
interaction, without destroying any symmetry. It can be considered as a natural infrared 
regulator. The adiabatic limit g ^ 1 must be performed at the end of the calculation in 
the right quantities where this limit exists. The existence of the adiabatic limit becomes 
a problem if the theory contains massless fields. 

The n-point operator-valued distributions r„ are the basic objects of the theory. They 
can be constructed inductively from Ti through a number of physical requirements, the 
most essential one being causality. Unitarity plays no essential role. Let the operator- 
valued distributions T„ be defined by 

°° 1 r 

S{g)-'^ = 1 + J2-, d^xi...d^XnTn{xi,...Xn)g{xi)...g{xn). (2.2) 

n=l 

Then one defines, for arbitrary sets of points X, Y in Minkowski space, the following 
distributions 

A!^{xi,...Xn) = Y.TnAX)Tn-nAy,Xn) (2-3) 

P2 

B!^{xi, . . . Xn) = {Y,Xn)Tn,{X), (2.4) 

P2 
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where the sums run over all partitions 

P2: {xi,...x„_i} = xuy, 

into disjoint subsets with |X| = ni, |y| < n — 2, and \X\ means the number of points in 
the set X. We also introduce 

Dn{xi,...Xn)=R'n-K- (2-5) 

If the sums are extended over all partitions P2, including the empty set X = 0, we obtain 
the distributions 

^n(^l) ■ ■ ■ = ^ ] Tn\ (-^)^n— ni {Yi ^n) = 

= < + T„(xi,...x„), (2.6) 

-^n(^l) ■ ■ ■ Xji) = ^ ] Tji—ji^ (Y, Xn)Tji^ (^) ~ 

= < + T„(xi,...x„). (2.7) 

These two distributions are not known by the induction assumption because they contain 
the unknown r„. Only the difference 

Dn = R'r,-A'^ = Rn- (2-8) 

is known. 

One can determine i?„ or An separately by investigating the support properties of the 
various distributions, this is the point where the causal structure becomes important. It 
turns out that Rn is a retarded and An an advanced distribution 

SUppEn C f+_i(a;n), SWppAn ^T~_i{Xn) (2.10) 

with 

^n-ii^) = • • • ^n-i) I Xj G V^{x),yj = 1, ... n - 1} (2.11) 

V^{x) = {y\{y- xf > 0, ±(yO - a;0) > 0}. (2.12) 

Hence, by splitting of the causal distribution (2.8) one gets i?„ (and A„), and r„ then 
follows from (2.7) (or (2.6)). The T„'s are well-defined time-ordered products. To carry out 
the splitting process, we write (2.8) in normally ordered form and then split the numerical 
distributions dn{x), where x = (xi — a;„, — x„). The causal splitting of d{x) can 

directly be done in momentum space by means of the following dispersion formula 

+00 ~ 

r{p) = ±- f dt ^, (2.13) 

2n J (i^iO)'^+i(l -t±iO)' ^ ^ 

—00 

which holds true for p G r+ (upper signs) or p G (lower signs) . The result for arbitrary 

p is obtained by analytic continuation, based on the fact that the retarded distribution 
f(p) is the boundary value of an analytic function, regular in R^""'^ -|- iT~^. This so-called 
central splitting solution is very convenient because it is obviously Lorentz covariant and 
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does not destroy any symmetry of the theory, uj is the singular order of d [2], in case of 
trivial splitting one has a; = — 1. 

The formulae (2.3-4) contain what is usually called Cutkosky rules. But we would like 
to emphasize that the sums run over all non-trivial partitions X UY. If X and Y are 
connected the term corresponds to an ordinary cut through the diagram. Otherwise the 
diagram is decomposed into more than two pieces. In this situation anomalous thresholds 
appear. But this only occurs in diagrams with more than three external legs. 

To calculate scalar diagrams with arbitrary masses, we start the inductive process from 
the following first order 

Ti(,t) = i : ip+{x)ip2{x) : A{x) - h.c. = -fi(x). (2.14) 

Here ipi, ip2 are charged scalar fields with masses m,i and m2, respectively, and A{x) is 
a neutral (self-adjoint) scalar field of mass ms. All fields are free fields satisfying the 
following commutation relations 

[^j{x),ip+{y)]=-iDm,{x-y), j = l,2 (2.15) 

[^'f\x),cpf+\y)] = -iD+.{x - y), (2.16) 

where Djn is the Jordan-Pauli distribution of mass m and (±) refers to the positive and 
negative frequency parts of the various quantities. A{x) fulfills the same commutation 
relations without the hermitian adjoint The commutator (2.16) gives the contraction 
in Wick's theorem, its Fourier transform is equal to 

DM = ^&{p')S{p'-m'). (2.17) 

The last vertex Xn plays a special role in the above equations (2.3-10). It is the splitting 
vertex, defining the edge of the causal cone. By translation invariance the numerical 
distributions only depend on the relative coordinates 

yj = Xj-Xn, j = l,...n-l. (2.18) 

The Fourier transform is always understood with respect to these relative coordinates 

d{p) = (27r)-2-+2 J d{y)e'^y d^i ■ • • d^yn-i- (2.19) 

Until now all n-th order distributions depend on n or n — 1 variables, the inner vertices 
are not integrated out. If the adiabatic limit exists, we can integrate the inner coordinates 
with g(x) = 1. In p-space this means that the inner momenta are put equal to 0. Then 
many terms in (2.3-4) vanish: 

Lemma 1. In the adiabatic limit only those partitions X L)Y contribute to (2-3) 
where X and {Y,Xn} contain external vertices, and similarly for R'^^ (2-4). 

Proof. Consider a partition where X contains no external vertex. Performing the 
contractions between X and {Y, Xn} with D+-distributions and transforming into p-space 
(2.17), we get a product of 9-functions 
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where all momenta add up to 0: 



Such a product is zero. 

The exists one serious problem, however. The central splitting solution (2.13) is only 
true if all momenta pj are inside the light cone. Therefore, strictly speaking, we cannot 
put the inner momenta equal to 0. But if the £)^-distribution is massive m > 0, the 
vanishing of the contribution of a wrong partition takes place for small enough inner 
momenta pj in y+, already. Then we can use (2.13) and take the limit pj 0. For this 
reason we always calculate with massive fields first. If the limit m — is required, it must 
be carefully performed by taking cancellations of infrared divergences between different 
terms into account. 

From lemma 1 it is clear that in diagrams with two and three external legs, the non- 
vanishing terms (in the adiabatic limit) correspond to ordinary cuts through the diagram. 
But in a four-legs diagram a pair of opposite legs can be in X and the other pair of opposite 
legs in {Y,Xn}- Then this decomposition is no longer a simple cut and an "anomalous 
threshold" appears. To our knowledge such a diagram was never computed by the naive 
"dispersive method", but in the causal theory this is no problem. The following second 
lemma further simplifies the later calculations. It is a consequence of parity- and time- 
reversal invar iance. 

Lemma 2. In a PT-invariant theory the numerical distributions d^{x) in Dn (2.5) are 
essentially real (i.e. up to an overall factor i) in momentum space: d^ip)* = d^{p)- 

Proof. We know that the causal ^-distributions are PT-invariant ([2], p.281). This 
implies for the numerical distributions: (i^(— x) = d^{x)*. The complex conjugate comes 
from the antiunitarity of time-reversal. After Fourier transformation this gives the desired 
result in momentum space. The overall factor i depends on whether the number of internal 
lines in the diagram is even or odd. 

We now come to the calculation of the third order vertex diagram (Fig. 2). To simplify 
the notation, we write the arguments in x-space without the dummy variable x. From 
(2.4) we have 

R's = T2{1, 3)fi(2) + T2{2, 3)fi(l) + Ti (3)^2(1, 2). (2.20) 
The first term herein contains a Compton subgraph 

R's, = i: ^+(1)L>^,(1 - 3)^2(3) :: ^(1)^(3) : {-i)A{2) : ^tM^) (2-21) 

where 

hF(^\ 

p2 _ ^2 _|_ 



D^ip) = ^':r.'.n (2-22) 



is the Feynman propagator. The product (2.21) is computed by Wicks theorem, restricting 
ourselves to those contractions which generate the vertex diagram: 

^31 = - : 9^^(lXi(l - 3)i?+,(3 - 2)D+3(1 - 2)991(2) : A(3). (2.23) 

Similarly we get for the other two terms in (2.20) 

R32 = --- <ft{^)D^,{3 - l)D^,{S - 2)L>+3(2 - 1)^1 (2) : ^(3), (2.24) 
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R's3 = -: vti^)Di, (3 - l)^+,(3 - 2)D^^ (1 - 2)^i(2) : A(3). (2.25) 
The anti-Feynman propagator D^^ is the complex conjugate of . It appears because 
we have used unitarity 72(1,2) = T2(l,2)"'". This is the only minor role which unitarity 
plays here. 

The result for A'^ (2.3) is obtained in the same way. Collecting the terms with field 
operators : (^2^(l)(^i(2) : A{3) in (2.5), the corresponding numerical distribution is 
given by 

d3(l,2,3) = -D^^il - 3)D+^{3 - 2)D+^{1 - 2) - D-^D^^D-^ 

AF 



+D~ + D~ D" + D'^ - D'^ D~ D 

' mi 7712 ms ' 7711 7712 7713 7711 7712 7713 7711 7712 7713 



(2.26) 

where we have used D~{x) = —D^[—x) and the arguments in the 6 terms agree with the 
first term. The 6 terms come from three cuts through the vertex diagram, not only one. 
To get contact with the convention in [2], we shall use the relative coordinates 

yi=xi- xa, 2/2 = 2:3 - X2, (2.27) 

and calculate the Fourier transform 



4(P,5) = (27r)-^y d3(yi,2/2)e^P^i+^^^^d^yiA2. 



(2.28) 



Then we arrive at 



d^{p,q) = {2^)- 



dk 



D-,^{p-k)Dt^{q-k)D^^^{k) 



-<i {P - k)Di, {q - k)D+^ (k) - D-^ {p - k)D^^ {q - k)D^^ {k)\ - [p ^ g, mi ^ ms] . 

(2.29) 

Up to the arbitrary masses and the imaginary parts, this agrees precisely with the result 
in QED ([2], eq.(3.8.28)). By the same techniques as in the QED case, we then find 

■ sgn Po , 



dzip, q) 



TT 



4(27r)f 



-e(P"-(mi + m2)")logi 



sgn go 



e(g2 _ + m^f) log2 +^^e(p2 - (mi + ma)^) loggj, (2.30) 



where P = p — q, N = (pq)^ — p^q^ and 

p2 + rnl-ml- pP(l + !I^)+^^f^ 



logi = log 



p!2 T I p!2 J 



pT- 



log2 = log 



logs = log 



p'^ — mf + ml 



^3 ^2 



47772 



— + ml — pq^l 



9 9 

-,2_^2 



NJ{1 



la\2 



Ami 



1 111 { 1 771?— 771, \ , rrf //-, 771?— 771, NO 4771, 

q^-mi + mi- pq[l - )+VNJ (1 ^)2 - 



— m2 + mg — pqll 



Ami 



(2.31) 



(2.32) 



(2.33) 



Because of lemma 2 we only need the real parts of the logarithms. The splitting of (2.30) 
by means of the central solution (2.13) is done later in (3.10). 
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3 Vacuum Polarization in Fourth Order 



Now we envisage the calculation of the diagram shown in Fig. la, which contributes to 
vacuum polarization in fourth order. We restrict ourselves to the case with vanishing 
'photon' mass ma and mass m of the 'electron'. Prom (2.4) we have 

E!^ = r3(l, 2, 4)Ti (3) + r3(i, 3, 4)ri (2) + r3(2, 3, 4)ri (1) 

+r2(l, 4)72(2, 3) + r2(2, 4)12(1, 3) + T2(3, 4)T2(1, 2) + ri(4)T3(l, 2, 3). (3.1) 

Note that we only consider terms with field operator : ^(2)yl(4) :. According to lemma 1, 
the first, the third and the fifth term in (3.1) vanish in the adiabatic limit. Furthermore, 
the second term gives the same contribution as R'^-j, and the same holds true for R'^ 
and i?4g. 

In x-space, the three-particle contribution Rl^ (Fig. 3a) is given by 

^^44 = ^44(1,2,3,4) :^(2)^(4) :, (3.2) 

r^4 = 2iDi^{2 - 3)D^{1 - 4)D+(1 - 2)D+(4 - 3)D+^{1 - 3) (3.3) 
The Fourier transform of (3.3) is (?/j = 

= 2(27r)-ii / dqdp' , 1 e{q^)6{q' - m') 

J [p2 + q) — — 7.0 [pi — q — p ) — + 

©(-P2 - - 9° - y°)'5((P2 +P3 + q + P' f - m2)e(p'°)5(p'' - ml) (3.4) 
which becomes in the adiabatic limit pi , Pa — ^ , p2 = P 

/\ \ 
dqdp' 2 ^ 7 7^2 JTT^ 
\P + q') — m"' — zO {q — p — p ) — + 

e{qo)S{q^ - m^)Q{p'o - qo)Si{p' - qf - m^)@{-po - p'^)5{{p+p'f - ml). (3.5) 

In fact, the calculation of (3.5) has already been performed in [2] by G.Kallen and 
A.Sabry. Our calculations confirm their result exactly, namely: 

"■^4^^^ = 4(^^®^~^°^®^^' ~ ^""'^ ^^-^^ 
B{z) = 3Li2(2;) + 2Li2(-z) + log z log(l - z) + log z log(l + z) + ^ log^ 2; - ^ + log z log 

(37) 

up to terms that vanish for ma 0. Here we have introduced the well-known dilogarithm 
Li2, and the variable 

z = (3.8) 
which varies from to 1 for p^ G [4m^, 00). 



8 



The two-particle contribution (Fig-3b) is 

rM = 2i(27r)-2 J dp' D+{-p')A^{p',p' - p)D+{p' - p). (3.9) 

Here, A3 is the vertex function. Since its first argument Ues in the backward hght-cone 
and the second one in the forward hght-cone, we are in the case of lower signs in (2.13), 
because of our convention (2.27). Then the retarded vertex function N!^^^ is given by the 
following dispersion integral: 

+00 
—00 

Again, the computation of (the real part of) can be carried out in a straightforward 
manner. As we know from (2.30), the vertex function consists of three different parts. 
The first one leads to an infrared finite contribution to 

1 TT^ 

Ci{z) = 2Li2(z) + 2 log z log(l -z)-- log2 z + —, (3.12) 

whereas the second and third term arc equal and contain an infrared divergent term 
~ log(m3/m), which cancels the infrared divergence in (3.7): 



r'l2{p) = J^^®i-Po)^(p' - 4"^') C2{z), (3.11) 



C2{z) = -U2{z) - logzlog(l -z) + ] log2 2 + ^ - logzlog ^. (3.12) 

4 6 m 

Now ^4 = r4 — 04 can immediately be written down if we note that 

= «4^(-p), i = i...,7. (3.13) 

4 The Splitting of d4(p) 

After having calculated the causal distribution 



diip) = - ^^^"'^^9 ^sgnpoQ(p^ - 4m2) J(z), 



J{z) = 4Li2(z) + 2Li2(-z) + 21og(z) log(l -z)+ log(z) log(l + z), (4.1) 

we must decompose it into retarded and advanced parts. The retarded distribution is 
given in the forward light-cone according to (2.13) 

+00 

dA{tp) 



nip) = ^ j dt- 
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— oo 

It is very convenient to introduce the new integration variable 

_ y/ty - ^f^p^ - 4m^ 1 x + l ,^2/2 tit _ 1 a; - 1 

This leads to (js G y+) 

1 

r4(p) = — ^— / dx{ h \ r-\ J{x)- (4.4) 

2(2Try0p^ J \ x + 1 x-z x-l zi ^ ' ^ ' 



The imaginary part zO in the denominator of (4.2) is included in p^ as discussed at the 
end of this section. 

Prom now on, r/^{j)) is considered as a function of z (3.8) 

. = i f6 - 2 - hJ^h\ = ^-^IZ^^_ (4.5) 
2\ y I ) 1 + ^1 - 4.m'/p^ ^ ' 

Then we note that z G [0, 1] for p^ G [4m^, cxd), and 

1 

R{z) = I dx{ ^ + + -^}j{x) (4.6) 

^ ^ J y x+\ x- \ z X- zi ^ ' ^ ' 



has the property R{z) = R{l/z). The integral 

1 

Ri = -2 J ^^J{x) = (4.7) 



is just a constant and can be calculated with the formulae given in [10,11]. The calculation 
of 



is most easily performed by first calculating the derivatives of R3 and R4: 







R',{z) = ^+ dx-P^. (4.9) 
z — I J 1 X — z 


This is a helpful trick, but all integrals coming up in (4.6) are also discussed in the 
literature mentioned above. We give the separate results for the two-particle and three- 
particle expressions: 

1 

/ ^^1 — TT + — \r + ^}^(^) 

J I- x + 1 x — \/z x — zJ 
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= 5Li3(2;) + 3Li3(-2;) - 3U{z) log z - 2U{-z) log z - ^ log' z log(l - z) 

1 1 tt' tt' 3 tt' 

-- log' z log(l + ^) - ^ log^ z + Y + -2) + Y log + -C(3) + Y log 2, (4.10) 

and 

1 

f dx{ ^ + + -^}c{x) 

J ^ x+1 x-l/z X- Zi ^ ' 


= Li3(^) - Li(^) log 5; - ^ log' z log(l - ^) + ^ log^ ^ 

+7r'log(l - z) - !^logz + ^C(3) - Y log 2, (4.11) 
where C {x) = Ci (x) + C2 (x) . Here we have introduced the trilogarithm, which is defined 

by 

U{x) 



Li3(z) = / dx=^. (4.12) 



The infrared divergent terms in B and C which cancel mutually in have already been 
omitted. Finally is given by 

''^'^ = -2(2;?k^^^^)' ^^-^'^ 



R{z) = GLhiz) + 3U3{-z) - AU2{z) log(-z) - 2Li2(-z) log(-z) 

- log'(-^) log(l -z)-^ log'(-z) log(l + z) + ^C(3). (4.14) 

If R{z) becomes complex the sign of the imaginary parts is determined by the iO in (4.2) as 
follows. For arbitrary time-like p, the iO goes over into iOpo (see (2.13) and [2], Chap. 3. 6). 
To obtain the full time-ordered distribution t^^p) wc have to subtract r'^{p)- This changes 
iOpo into iO again. Consequently, t4{p) is given by (4.13) with j?' substituted by p' -|- iO, 
i.e. z ^ z + iO. This fixes the signs of the imaginary parts at the logarithmic cuts in 
(4.14). For space- like p, t4{p) is simply given by (4.13) because r4(p) now vanishes and 
R{z) is real. 



5 The Electron Propagator in Fourth Order 

After having demonstrated the instructive example of vacTium polarization, we proceed 
now with the discussion of the electron propagator in fourth order (Fig. lb). Again, we start 
from the general expression (3.1) and consider only terms with field operator : (p(2)(p~^{4) :. 
Just as in the case of vacuum polarization, the first, the third and the fifth term in (3.1) 
vanish in the adiabatic limit, and the second and the seventh term 

R'42 = Ts{l,3,A)fi{2), i?l7 = Ti(4)r3(l,2,3) (5.1) 

give the same contribution to the propagator. But we have to distinguish the different 
cuts in i?44 and R'^q. 
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In x-space, the three-particle contribution in R'^q with two photons and one electron 
as intermediate state (Fig.4) is given by 

i?'l6 = ^'16(1,2,3,4) :(p(2)(^+ (4):, (5.2) 

r% = iDi{3 - 4)Di^{l - 2)Z)+(3 - l)D+^{4: - 1)0+^(3 - 2) (5.3) 
The Fourier transform is after the adiabatic limit 

1 1 



r%{p) = {2it)-^^ J dqdp'- 



e(go)<5(g' - mi)e(p'o - q^)5{{p' - qf - mi)e(-po - p'o)<5((P + p'f - rr?). (5.4) 
Straightforward calculation of (5.4) leads to the simple result 

""'^6 = 8(^^®^~^°^®^^' " m')B{x% (5.5) 

S(x2) = iLi2(x2) + ilog(x2)log(a;2-l)-^ , x'=p^/m\ (5.6) 

The two-particle contribution can be evaluated without any problem with the aid of (2.30). 
We therefore quote only the result: 

i?'42 = r'42(l,2,3,4) : ^(2)<p+(4) :, (5.7) 
r'42(p) = i(^^0(-^^o)e(p2 _ ^2^c{x\ (5.8) 

C{x^) = \u^{x^) - (5.9) 

In (5.6) and (5.9) only the real parts of the dilogarithms contribute. 

Since wc have to treat the thrcc-particlc contribution with three electrons as intermedi- 
ate state (r'44 and t'^q) separately, we already give the splitting results for the expressions 
obtained so far. The splitting procedure leads to the same type of integrals as in the case 
of vacuum polarization in Sect. 4, hence we refer again to [10,11]. We have for p in the 
forward light-cone and p^ > 



1 1 



di{p) = -^^^^^sgnpo0(p -m^)[3Li2(x^) + logx^log(x^ - 1) - yj, (5.10) 

+00 

nip) = f / dt-^lM^ (5.11) 

—00 

~ 8(27r)iV [^Li3(l-x^)-31og(l-x^)Li2(l-x^)-logx^log^(l-x^)-4C(3)]. (5.12) 
The calculation of the second diagram in Fig.4 

1 1 



r'Up) = (277)-'' J dqdp' 



{p + q)^ ~ "^3 ~ {q ~ P ~ P')^ ~ ™-3 + ^0 
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e{qo)5{q'' - m2)e(p'o - qo)S{{p' - qf - rr?)Q{-po - p'o)S{{p + p'f - rr?). (5.13) 

is difficult and requires extra consideration. r''\^{p) turns out to be infrared finite, so we 
may drop the small photon mass in the fofiowing. Making use of all ^-distributions, the 
integral 

1 = j dqe{qo)5{q'-m')Q{p',-qo)5{{p'-qf-W?) ^^^^^,_^^ (g - + ,0 ^^'^^^ 

can be transformed into 

^ ^ _.e(rt)e(p"-W) 7 1 

Jip'pY -p''^p^ + m''){x + pp) 

where 



p = -p-p' , xi,2 = ^p'p T [i^pf - Vy 1 - (5.16) 



Then we obtain 



x^,2 = l{y-m')T^^, ^={y'-p'm')Jl- , J"^^ , ■ (5.18) 
2 2 y + m"' — 2y 

Introducing x = (y + z — m^)/2, d2(p) = 2(r'44(^?) — a'%^(j))) becomes 

i{p2_3m2) +^ 
11 f f 1 

t^2(p) = - ..^ x9 — sgnpo 9(p^-9m^) / / dz— — — — ^ — ^ (5.19) 

4(27r)" ^ J {y + z + m'^jyy — z + m'^) 

At this stage, we will not proceed the same way as in the case of vacuum polarization. 
We first apply the splitting formula to (5.19) 

4iP) 



"4(27r)iV / ^(1^7+^ / / '^\y + z + m?){y-z + 'm?y ^^'^^^ 



?(5,y) = \/y'-p'm2sjl + - ^"^^ , , (5.21) 

2y — — p^s 
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and perform a partial integration with respect to s. This leads to the following integral: 

+00 



i(p2s-3m2) 



r2{p) 



2(27r)iV 



ds log(— 



s-iO 



dy 



A 

ds 



(5.22) 



+00 J/2 ^1 

2 2/2 -'-/2ii\ ^/2 o\ 

X =p /m , yi = -[X s + 1), y2 = -{xs-3) 



We obtain after some simple manipulations 



i y2-y 
yi-y' 



r2 



I 



2(27r)i0m2 J {x^s-l^ 

9/a;2 



J/2 

/ 



dy 



^{y^ - x^s){y - yi){y - ^2) 



{x^s-l)y-^{x^s + lf + 1 



(5.23) 



(5.24) 



(5.25) 



The last integral can be expressed by complete elliptic integrals of the first and third 
kind. It is 



2/2 

/ 



dy 



Viy' -x''s){yi-y){y2-y) ^'(y^ _ V^)(y2 + V^) 



where 



and 



, 2 _ (2/2 - V x'^s){yi + V x^s) 



2/2 

/ 



{yi - Vx^){y2 + Vx^s) 



ydy 

V'(y^ - x2s)(j/i - y)(y2 - y) 



K(A;), (5.26) 



(5.27) 



(yi - Vx^){y2 + WsFs) 



2Vx2^n(a^ fc) - \/x^sK{k) 



a 



,2 _ 2/2 - VX^S 



(5.28) 



(5.29) 



2/2 + va;^s 

Introducing A = Vx^s G [3, 00), we find that the modulus k and the parameter are 
related by the following identities: 



(A-3)(A + 1) 2 (A-3)(A + 1)^ 

a = TTT-: — , K — 



(A + 3)(A-1) 



(A + 3)(A-1)3- 



(5.30) 
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Then it is in fact possible to express 11(0^, k) by K{k): 

U{a\k) = ^^K{k). (5.31) 

This has already been observed by A.Sabry [12]. Our relation follows from his equation 
(85) by the substitution A — 1/A. Finally we arrive at a remarkably simple expression for 

oo I 

i f A3 /A I 2 s 

''^^^ = 6(2^ J (A - 1)^(A + 1) V A^ ^°g(Ti:7^)^(fe)- ^5-32) 

9/a;2 

So far we have calculated the retarded distribution for time-like = p^/rn^. By the 
same argument as given at the end of the last section, this also gives the time-ordered 
distribution t2{p) for arbitrary p^. To get the full fermion propagator, one must add ti{p) 
given by (5.12) with ^ x'^ + iO. The integral in (5.32), which can easily be computed 
numerically, is in agreement with eq. (30) of Broadhurst [4]. Our funny normalization 
factors are the correct ones for the calculation of S-matrix elements according to (2.1). 
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Figure Captions 



Fig.l The two- loop master diagrams: 

(a) the 'photon' propagator; 

(b) the 'electron' propagator 

Fig. 2 The vertex diagram with arbitrary masses 

Fig.3 (a) The three-particle cut, (b) the two-particle cut for fourth order 

vacuum polarization 
Fig.4 The two different three-particle cuts in the electron propagator 
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